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SECTION  1 


RESEARCH  OBJECTIVES 


1.1  OVERALL  OBJECTIVE 

The  objective  of  this  research  program  is  to  develop  an  understanding  of 
oxidizer  particle  size  effects  on  the  non-steady  combustion  behavior  of  compo¬ 
site  solid  propellants.  This  understanding  is  expressed  by  the  development  of 
suitable  analytical  models  to  describe  known  or  experimentally  observed  pheno¬ 
mena.  A  combination  of  theoretical  and  experimental  tasks  is  performed  toward 
the  accomplishment  of  this  objective.  Results  are  disseminated  by  technical 
presentations  and  publications,  and  through  a  variety  of  interchange  functions 
with  government  laboratories  and  contractors.  Demands  of  future  propulsion 
systems  and  the  viability  of  combustion  tailoring  through  control  of  oxidizer 
particle  size  distribution  make  it  imperative  that  the  effects  of  size  distri¬ 
bution  be  understood. 


1.2  SPECIFIC  OBJECTIVES  FOR  FISCAL  YEAR  19B2 

Particular  attention  has  been  given  to  the  oroblem  of  combustion  insta¬ 
bility,  which  is  promoted  by  the  devel ooment  of  improved  rocket  motor  capabili¬ 
ties  for  future  missile  systems.  The  specific  objectives  are  listed  as 
follows: 

(1)  Develop  an  analytical  model  for  the  velocity-coupled  combustion 
response  function,  accounting  for  effects  of  ammonium  perchlorate 
(AP)  composite  propellant  heterogeneity. 

(2)  Refine  and  apply  experimental  techniques  utilized  during  FY  ig81  to 
characterize  the  composite  propellant  heterogeneity. 

(3)  Begin  analytical  work  to  characterize  the  high  freguency  pressure- 
coupled  response  properties  of  minimum  smoke  (nitramine/energetic 
binder)  propellants. 

A  literature  survey  on  analytical  response  function  models  which  account 
for  effects  of  AP  particle  size  (i.e.  ,  heterogeneity)  was  completed  during 
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Ing  1981  for  the  pressure-coupl ed  response  function  (Refs.  2,  3).  This  model 
accounts  for  the  heterogeneity  in  two  ways:  affects  of  AP  particle  size  on 
ballistics  properties  which  enter  into  the  response  function  theory;  and  compo¬ 
sitional  fluctuations  of  frequencies  associated  with  characteristic  dimensions 
of  the  heterogeneity.  The  former  is  in  the  nature  of  classical  response  func¬ 
tion  theory;  the  latter  incorporates  a  new  mechanism  by  which  the  combustion 
process  can  drive  an  instability.  Open  publication  of  this  work  was  recently 
accomplished  {Ref.  4).  The  work  planned  for  FY  19B2  consisted  of  moving  on  to 
the  suhject  of  the  velocity-coupled  response.  The  presence  of  crossflow  and 
fluctuations  in  crossflow  velocity  are  important  aspects  of  combustion  instabi¬ 
lity,  and  more  difficult  to  treat  than  pressure  fluctuations  (Ref.  5).  Many 


believe  velocity-coupling  to  be  an  essential  ingredient  of  non-linear  instabi¬ 
lity,  by  which  an  instability  becomes  most  destructive.  As  a  first  step,  howe¬ 
ver,  the  analytical  work  would  remain  in  the  regime  of  linear  instability.  The 
linearized  problem  is  simpler  to  treat  within  the  current  state  of  knowledge, 
and  is  relevant  to  the  question  of  stahility  margin  (Ref.  6).  It  is  considered 
that  much  can  he  learned  from  the  properties  of  the  linear  velocity-coupled 
response.  It  is  planned  to  begin  work  on  the  non-linear  response  in  1983. 

Three  types  of  experiments  have  heen  explored  to  characterize  the  heter¬ 
ogeneity  of  composite  propellants  and  its  effects  on  combustion  behavior.  The 
first  experiment,  which  is  non-destructive,  is  to  perform  energy  dispersive 
analyses  of  x-rays  (EDAX)  scans  of  propellant  samples,  and  Fourier  decomposi¬ 
tions  uf  the  measured  fluctuations  in  chlorine  intensities  representative  of 
local  AP  concentrations.  The  second  experiment  is  to  use  the  JPl  microwave 
hurner  (Ref.  7)  to  measure  dynamic  burning  rates  of  the  same  propellants  at 
constant  pressure,  and  to  analyze  the  data  for  frequency  components.  The  third 
experiment  is  to  use  the  burner  to  measure  pressure-coupled  response  functions 
of  these  propellants  and  determine  peak  response  frequencies.  It  is  desired  to 
seek  the  consistent  appearance  of  any  preferred  frequencies,  and  to  relate  the 
observed  frequencies  to  the  heterogeneity  of  real  propellants.  Feasibility  of 
the  approach  was  established  in  FY  1981  work  (Refs.  2,  8),  but  technique  re¬ 
finements  are  needed  to  extract  the  preferred  frequencies,  and  more  experiments 
are  needed  to  establish  their  consistency.  Plans  for  FY  1982  included  continu¬ 
ation  of  this  experimental  work. 

Minimum  smoke  propellants  have  exhibited  a  propensity  toward  high 
frequency  tangential  mode  instability  (Ref.  9).  The  high  frequency  regime 
introduces  the  question  of  gas  phase  time  lag  as  an  important  contributing 
mechanism,  and  therefore  the  need  to  relax  the  guasi-steady  gas  assumption  com¬ 
monly  used  in  response  function  theories.  One  work,  by  Tien,  has  relaxed  this 
assumption  (Ref.  10),  but  the  results  are  suspect  because  of  the  prediction  of 
a  negative  response  over  a  bread  freguency  range  that  is  contrary  to  experi¬ 
ence.  Another  interesting  aspect  of  minimum  smoke  propellants  is  that  their 
combustion  behavior  exhibits  the  characteristics  of  homogeneous  propellants 
(Ref.  11).  Therefore,  an  approach  that  neglects  heterogeneity  but  relaxes  the 
quasi-steady  gas  assumption  appears  viable.  It  was  planned  to  review  the  basis 
of  Tien's  work  with  intentions  to  revise  that  work  as  necessary  and  study  the 
resulting  properties.  The  pressure-coupled  response  properties  of  minimum 
smoke  propellants  are  important  because  high  frequency  stability  additives  may 
be  precluded  or  ineffective  (Ref.  9). 
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SECTION  2 


STATUS  OF  THE  RESEARCH  EFFORT 


2.1  VELOCITY-COUPLER  RESPONSE  MOOELING 

2.1.1  Analysis  of  the  Linear  Vel ncity-Coupled  Response  Function 

A  theoretical  analysis  of  the  linear  velocity-couoled  response  function 
has  been  completed.  This  work  coniplements  the  previnus  analytical  developments 
for  the  pressure-coupled  response  functinn,  including  a  compositional  response 
mechanism  to  account  for  composite  ornpellant  heterngeneity  (Refs.  2-4).  There 
are  two  aspects  of  the  prohlem  of  linear  velocity-coupl inq.  First,  as  pointed 
out  hy  Lengelle  (Ref.  12),  is  the  effect  of  mean  crossflnw  velocity  on  the 
various  response  function  parameters.  Second  is  the  combustion  response  to 
fluctuations  in  the  velocity,  which  is  the  velocity-coupled  response  function. 
Roth  aspects  are  of  interest  in  this  study. 

The  basic  mechanisms  are  time  lags  in  the  solid  phase  driven  by  pressure 
and  velocity  perturbations  in  the  qas  phase.  Added  to  these  perturbations  are 
compositional  fluctuations  which  originate  in  the  solid  phase;  this  is  the 
effect  of  heterogeneity.  The  analysis  is  limited  to  conditions  under  which 
quasi-steady  gas  assumptions  and  principles  of  linearization  apply. 

The  effect  of  crossflow  velocity  basically  involves  a  boundary  layer 
heat  transfer  mechanism,  by  which  the  static  heat  feedback  from  the  combustion 
zone  is  modified  and  augmented.  Under  steady-state  conditions,  the  phenomenon 
is  comnonly  referred  to  as  "erosive  burning."  Where  the  flow  has  an  oscillat¬ 
ing  component,  as  during  unstable  burning,  the  combustion  response  to  this 
oscillating  component  is  sometimes  referred  to  as  the  "erosive  response,”  but 
is  more  often  called  the  'velocity-coupled  response.”  Velocity-coupling  itself 
may  involve  more  than  combustion  response  (Ref.  5),  hut  this  research  is  limit¬ 
ed  to  the  question  of  combustion  response.  A  numher  of  theoretical  analyses 
describing  the  houndary  layer  heat  transfer  mechanism  in  the  context  of  snlid 
propellant  combustion  have  been  puhlished  in  recent  years  (Refs.  12-21).  There 
is  much  disagreement  in  describing  the  mechanism,  and  the  problem  has  not  been 
resol ved. 

The  approach  adopted  in  this  wnrk  has  been  to  consider  the  boundary 
layer  interaction  mechanism  as  a  separate  prohlem,  and  concentrate  on  the  com¬ 
bustion  response  to  whatever  the  perturbing  function  is  that  is  transmitted. 
The  boundary  layer  interaction  mechanism  is  a  suhject  of  continuing  work  by 
others  (Refs.  22,  23)  and,  when  sufficient  progress  is  achieved,  the  intention 
is  to  combine  the  elements  into  a  fully  coupled  model  that  would  ultimately 
entail  nonlinear  effects  (Refs.  24,  25)  and  probably  the  core  gas  dynamics 
(Ref.  23)  as  wel 1 . 

In  view  of  the  above-mentioned  uncertainty,  and  the  desire  to  concen¬ 
trate  on  the  combustion  response  (with  emphasis  on  AP  particle  size  effects), 
it  was  assumed  that  influences  of  crossflow  velocity  could  be  represented  by  an 
addition  to  the  surface  energy  balance  taking  the  form  of  Saderholms’  law  (Ref. 
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27).  This  particular  law  was  selected  because  it  has  a  qood  history  of  corre- 
latino  erosive  burninq  data  over  a  wide  ranqe  of  conditions  (Refs.  28,  29), 
appears  to  have  the  orooer  interplav  between  controlling  variables  that  any 
mechanistic-based  theory  would  have  to  match,  and  has  been  olven  mechanistic 
support  in  the  work  of  Kuo,  et  al.  (Ref.  30).  The  Saderholm  expression  would 
be  of  the  form 


E  =  Kj  (u0/ut) 


\o)N’ 


(I) 


Assuming  the  erosive  threshold,  ut ,  to  increase  with  increasing  burninq  rate 
and  decreasing  particle  size  (Ref.  14),  Eq.  (1)  becomes  the  adoitive  enerqy 
transport  term 


Nu  Np  N0  Nr 

E  =  (u0)  (P0)  (D0)  /  (ro) 


(2) 


This  equation  was  optimized  to  fit  the  Cohen-Strand  model  of  steady-state  burn¬ 
ing  (Ref.  31)  to  all  of  the  experimental  data  of  Kinq  (Ref.  18).  The  resulting 
standard  deviation  between  calculated  and  experimental  burninq  rates,  r,  when 
compared  as 


x  100 


was  24X.  This  is  considered  quite  qood  when  compared  on  this  basis  (Ref.  32), 
because  it  is  comparinq  the  incremental  burninq  rates  rather  than  the  burninq 
rates  themselves.  Also,  it  better  reflects  the  error  between  the  calculated  E 
and  that  value  of  E  required  in  the  model  to  match  the  experimental  result 
exactly. 


It  should  be  restated  that  the  significance  of  E  in  this  approach  is  as 
part  of  a  boundary  condition,  to  represent  a  transport  mechanism  that  has  not 
yet  been  conclusively  quantified  in  independent  theoretical  analysis.  As  such, 
its  semi -empi n cal  origin  does  not  render  the  approach  unscientific  or  without 
physical  significance.  For  example,  it  is  of  interest  that  the  value  of  the 
particle  size  -  dependence,  Nq,  is  a  number  close  to  zero  ( Ng  =  0.15).  This 
means  that  the  effect  of  a  coarse  particle  is  largely  due  to  the  fact  that  the 
blowing  (r0)  tends  to  bp  low,  and  not  due  to  the  fact  that  there  is  a  big 
particle  nicking  up  into  the  boundary  layer.  Also,  the  value  of  the  velocity- 
dependence,  Nu,  is  a  number  close  to  0.8  (Nu  =  0.87).  Such  a  dependence  is 
mindful  of  classical  convecti ve-heatinq. 

With  E  thus  established  *..•  added  to  the  model,  a  perturbation  analysis 
provides  the  velocity-coupled  response.  The  three  response  function  elements 
which  have  been  derived  in  the  course  of  this  work  may  be  expressed  as: 


2-2 


{3} 

(4) 


Rp/np  =  1/1 

Ra/n a  3  [1  -  (V na)  °p  <TW  -  T0)  (1  -  1/X)]/L 


^u^nu  "  (Vuo^L 


(3) 


where 


L  *  1  -  op  (T*0  -  T0)  (1  -  1/X)  *  aj  (X  -  1) 


°T  * 


n 


n 


n 


P,  a,  u 


P,  a,  u 

a.  u,  T0 

P.  u.  T0 

P.  T0 


(6) 

(7) 

(«) 

(9) 

(10) 


(11) 


The  real  parts  of  Eqs.  (3)  and  (4)  combine  to  achieve  the  overall  driv¬ 
ing  pressure-coupled  response:  the  Imaginary  parts  of  Eqs.  (4)  and  (5)  combine 
to  achieve  the  overall  driving  velocity-coupled  response. 


2.1.2  Parametric  Studies 

Series  of  computations  were  per'ormed  over  ranges  of  the  controlling 
comb- *.t Ion  parameters,  as  reported  previously  (Refs.  2-4),  but  now  the  mean 
flow  velocity  *s  added  as  a  parameter  and  the  velocity-coupled  response  para¬ 
meters  are  added  as  calculated  dependent  variables.  Calculations  were  made  at 
velocities  of  200,  400  and  600  m/sec.  Results  at  0  m/sec,  where  the  velocity- 
coupled  response  would  be  0  (u^>0),  were  reported  previously.  The  review  of 
the  results  may  be  highlighted  as  follows. 
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With  respect  to  the  influence  of  mean  flow  velocity  on  the  controlling 
ballistics  properties.  Increasing  u0  tends  to 

-  increase  burning  rate,  r0 

-  Increase  pressure  exponent,  np 

-  reduce  burn  rate-temperature  sensitivity,  op 

-  increase  surface  temperature,  (Tw  -T0) 

o 

-  decrease  surface  temperature-temperature  sensitivity,  oy,  or 
cause  It  to  go  through  a  minimum  at  an  intermediate  velocity 

-  decrease  concentration  exponent,  na 

-  decrease  velocity  exponent,  nu 

W<th  respect  to  the  influence  of  mean  flow  velocity  on  the  pressure- 
coupled  ,,»ponse  function  elements,  Increasing  u0  tends  to 

-  change  the  peak  value  of  Rp  in  a  non-systemati c  way,  but  the 
changes  are  generally  small;  the  larger  peak  values  at  u0  =  0 
tend  to  decrease  with  increasing  u0,  and  this  change  is  more 
substantl al 

-  decrease  the  peak  value  of  Ra,  the  largest  decreases  occurring 
for  the  largest  peak  values  at  u0  =  0 

-  Increase  peak  response  frequencies,  generally  due  to  the  increase 
in  r0’,  but  sometimes  there  also  Is  a  step-up  in  che  peak 
dimensionless  frequency  i2;  usually  the  peak-response  12  decreases 
very  sliyhtiy  with  increasing  uQ. 

These  results  are,  of  course,  colored  by  the  particular  combustion  model 
that  is  used,  but  they  do  appear  reasonable.  Combining  the  pressure  -coupled 
response  elements,  increasing  the  crossflow  velocity  will  tend  to  be  stabiliz¬ 
ing  at  lower  frequencies  (near  and  below  the  peak  response  frequencies  at 
u0  3  0),  and  to  shift  the  response  function  curve  to  higher  frequencies, 
where  It  will  eventually  be  destabilizing.  Effects  of  A P  particle  size  are  as 
reported  previously  (Refs.  2-4). 

Results  for  the  properties  of  the  velocity-coupled  response  are  fascin¬ 
ating.  First,  the  typical  appearance  of  a  velocity  response  function  curve  is 
displayed  in  Figure  1.  The  Imaginary  part  of  the  velocity-coupled  response  Is 
the  component  that  Is  pertinent  to  the  combustion  instability  driving.  There 
are  two  regions  of  interest:  at  low  frequency,  burn  rate  fluctuations  are 
phase-leading  and  the  Imaginary  part  is  positive;  at  high  frequency,  burn  rate 
fluctuations  are  phase-lagging  and  the  imaginary  part  is  negative.  This  dis¬ 
tinction  is  very  Important  to  the  grain  design  of  a  rocket  motor  because  it 
determines  whether  velocity-coupling  adds  or  subtracts  acoustic  energy  at  par¬ 
ticular  locations  (Refs.  33,  34),  Three  key  frequencies  are  noted  on  the  plot; 
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that  corresponding  to  the  positive  peak  response:  (ftp),  that  corresponding  to 
the  negative  peak  response  (ft^),  and  that  corresponding  to  the  sign  change 
(ST*).  It  is  a  typical  result  that  the  negative  values  reach  somewhat  larger 
absolute  magnitudes  than  the  positive  values,  and  that  the  negative  part  exhi¬ 
bits  a  relatively  broad  and  flat  peak  region  compared  to  the  positive  part. 

The  Imaginary  part  of  the  comp- si tlonal  response,  Ra^»  exhibits 
qualitatively  similar  features  and  properties  as  the  imaginary  part  of  the 
velocity  response.  This  was  also  true  in  comparing  ths  real  parts  of  the 
compositional  and  pressure  responses,  as  reported  previously  (Refs.  2-4).  How¬ 
ever,  whereas  the  compositional  response  always  achieved  magnitudes  far  larger 
than  the  pressure  response,  such  that  it  could  be  inferred  that  the  heterogene¬ 
ity  would  dominate  pressure-coupled  driving,  the  magnitudes  of  the  composition¬ 
al  and  velocity  responses  are  more  comparable.  The  magnitudes  of  the  composi¬ 
tional  response  tend  to  exceed  the  velocity  response  at  lower  AP  concentra¬ 
tions,  lower  pressures,  and  wit!-,  finer  HP  particle  sizes.  The  magnitudes  of 
the  velocity  response  tend  to  exceed  the  compositional  response  at  higher  W 
concentrations,  higher  pressures  and  with  coarser  AP  particle  sizes.  There¬ 
fore,  with  regard  to  velocity-coupling,  the  heterogeneity  is  important,  though 
not  as  dominating  an  Influence  as  with  pressure-coupling. 

Another  consequence  of  relatively  large  values  of  both  the  compositional 
and  velocity  responses  is  that  the  velocity-coupled  response  can  tend  to  have 
larger  values  than  the  pressure-coupled  response.  Experimental  measurements  or 
deductions  of  absolute  values  in  excess  ef  1C  have  been  frequent  and  often 
viewed  with  suspicion.  This  *urx  supports  such  large  values,  and  does  so  with 
confidence  because  theoretical  and  experimental  values  of  the  exponents  (for 
pressure,  np,  for  composition,  na,  and  for  velocity,  n„)  are  in  fairly  good 
agreement  and  because  values  of  nu  are  generally  much  larger  than  values  of 

np  • 


Qualitative  effects  of  the  temperature-sensitivity  parameters,  op  and 
of,  upon  the  velocity  response  are  the  same  as  upon  the  pressure  response. 
Increasing  op  tends  to  increase  the  peak  values  arc  narrow  or  sharpen  the 
peak  regions.  Increasing  oj  tends  to  decrease  the  peak  values  and  narrow  the 
peak  regions.  The  response  is  proportional  to  the  exponent.  Effects  of  AP 
particle  size  upon  nu  are  not  systematic,  but  nu  generally  tends  to  in¬ 
crease  with  increasing  AP  content,  increasing  pressure  and  decreasing  velocity. 

Figure  2  shows  the  effects  of  AP  particle  size  and  pressure  upon  the 

nhaco-l  parll  nn  no»V  values  r>f  H10  imaninarv  nart  of  the  velocity  resnnnse. 

I  •  —  —  —  u  r  ™  ■  u  1  ■  ■  v  * 

These  calculations  were  made  for  88?  AP  and  a  mean  flow  velocity  of  200  m/sec, 
but  the  results  typify  the  qualitative  behavior.  The  most  striking  feature  is 
the  sharp  increase  with  pressure  in  the  coarse  AP  regime,  and  the  very  large 
peak  values  attainable.  In  general,  the  peak  response  values  Increase  with 
pressure  and  particle  size  although  an  intermediate  maximum  appears  at  an 
intermediate  particle  size.  Figure  3  shows  the  effects  of  AP  particle  size  and 
crossflow  velocity  at  constant  pressure.  The  striking  feature  here  is  the 
sharp  increase  with  decreasing  velocity  in  the  coarse  AP  regime.  In  general, 
the  peak  response  values  increase  with  decreasing  velocity  over  the  range  of 
velocities  considered.  Effects  upon  the  phase  lagging  peak,  and  upon  the  ana¬ 
logous  peaks  for  Ra'1)  are  qualitatively  the  same.  More  calculations 
with  further  decreases  in  crossflow  velocity  would  be  ot  interest. 
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PARTICLE  SIZE,  MICROMETERS 


The  three  key  frequencies  tend  to  increase  with  variables  which  bring 
about  an  increase  in  burning  rate:  increased  pressure,  velocity  and  AP  con¬ 
tent,  and  decreased  particle  size.  An  interesting  result  is  a  step-up  in  the 
values  of  U  that  sometimes  occurs  in  the  course  of  each  of  these  variations  in 
the  independent  variables.  Combined  with  the  normal  burn  rate  increase,  this 
step-up  in  the  SI  causes  a  relatively  large  and  sudden  shift  in  the  frequen- 
c*e^.  The  reason  for  this  hehavior  is  not  clear,  and  merits  further  study 
because  the  frequency  separating  phase  lead  and  phase  lag  response,  in  particu¬ 
lar,  is  of  special  importance  and  whatever  causes  it  to  change  sharply  would  be 
of  interest. 

Perhaps  the  most  important  highlight  is  that  the  velocity-coupled  res¬ 
ponse  is  largest  at  combinations  of  coarse  AP/low  burn  rate,  high  pressure  and 
low  mean  flow  velocity  (but  in  excess  of  threshold  velocity).  This  is  very 
gratifying  because  it  is  entirely  consistent  with  experience.  Figures  2  and  3 
would  suggest  that  400  ^m  AP  be  avoided,  especially  due  to  the  fact  that  this 
material  varies  from  lot-to-lot  and  often  contains  tails  of  very  coarse  parti¬ 
cles  in  the  size  distribution.  Other  figures  of  this  type  could  be  generated 
to  furnish  developmental  guidelines,  but  it  is  cautioned  that  specific  results 
are  suhject  to  uncertainties  in  the  comhustion  model  and  that  the  heterogenei¬ 
ties  of  propellants  must  be  determined  hefnre  this  theoretical  work  can  be  pro¬ 
perly  applied.  Nevertheless,  there  is  much  food  for  thought  and  future  work 
could  profit  from  delving  further  into  the  reasons  for  the  calculated  trends. 

2.1.3  Pulse-Triggered  Nonlinear  Instability:  Review  of  CARDE  Experiments 

A  classical  series  of  experiments  was  carried  out  by  Rrownlee  (Ref.  35) 
in  which  the  ability  to  trigger  axial  mode,  nonlinear  instability  by  pulse-in¬ 
jection  in  tubular-grain  rocket  motors  was  mapped  as  a  function  of  pressure  and 
geometric  variahles.  He  correlated  the  data  in  terms  of  a  critical  pressure- 
Kn  relationshio,  and  mapped  a  stability  boundary  as  Kn  versus  port  diame¬ 
ter.  These  relationships  were  found  to  be  dependent  upon  propellant  formula¬ 
tion  in  a  given  motor,  and  upon  the  type  of  motor  (Refs.  35-37).  Significant¬ 
ly.  slower  burning-rate  propellants  were  always  more  susceptible  to  this  type 
of  triggering  in  the  experiments  cunuuuted.  In  view  of  the  gpnpral ly-held 
belief  that  velocity-coupling  is  responsible  for  this  type  of  instability  (Ref. 
5),  a  review  of  the  work  was  deemed  appropriate  as  part  of  the  evolution  of 
this  research  into  the  area  of  nonlinear  velocity-coupling. 

An  important  characteristic  of  the  Brownlee  experiments  is  that  it  was 
conducted  by  reducing  the  nozzle  throat  area  on  successive  motors.  This  served 
to  increase  the  ranges  of  Kn  and  pressure  encountered  from  motor  to  motor, 
using  a  progressi ve-hurni ng  tubular  grain.  Relow  a  certain  pressure/Kn,  the 
instability  could  not  be  triggered.  Above  that  pre$sure/Kn,  the  triggered 
instability  progressively  increased  in  severity  with  increasing  pres$ure/Kn. 
By  reducing  the  throat  size  from  motor  to  motor,  the  instability  could  be  trig¬ 
gered  earlier  in  the  firing  and  stronger  instability  could  be  excited  later  in 
the  firing. 

The  Rrownlee  data  were  examined  to  see  if  the  triggering  could  be  asso¬ 
ciated  with  the  achievement  of  some  critical  low  axial  velocity  in  the  port. 
It  was  noticed  that  the  results  might  also  be  interpreted  on  the  hasis  of  a 
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critical  J,  because  the  smaller  the  throat  the  smaller  the  port  at  which  trig¬ 
gering  could  occur.  Sufficient  information  was  provided  to  enable  calculations 
of  the  port  exit  velocity  correspond! ng  to  each  triggering  boundary  point.  For 
ports  ranging  from  2.64-3.77  cm.  diameter,  Kn  varying  between  237  and  332, 
and  pressures  ranging  from  2.46-3.84  MPa,  at  triggering,  the  exit  velocities 
were  computed  to  vary  from  111-116  m/sec.  This  near-constant  value  is  consi¬ 
dered  to  be  significant,  although  the  significance  at  this  time  can  only  be 
speculated.  It  might  be  an  Indication  of  the  erosive  threshold  velocity.  It 
might  be  related  to  the  pulse  strength  In  terms  of  the  ability  of  the  compres¬ 
sion  wave  reflection  to  stop  or  reverse  the  flow,  or  cause  the  flow  to  fall 
below  the  erosive  threshold.  In  any  event,  this  observation  should  be  noted 
for  future  work. 

Later  experiments  (Refs.  36,  37)  were  not  so  extensive  as  to  permit  the 
evaluation  of  the  consistency  of  the  critical  velocity  for  a  given  propellant 
and  motor,  but  a  value  for  the  critical  velocity  could  be  determined.  Results 
for  five  propellants  are  tabulated  below,  relating  the  port  exit  velocity  to 
the  propellant  burn  rate  at  the  triggering  boundary  point. 


Burn  Rate  (cm/sec) 


Port  Exit  Velocity  (m/sec) 


0.31 

91 

0.40 

114 

0.63 

202 

0.70 

193 

0.82 

216 

With  one  exception,  there  is  a  good  relationship  hetween  the  burn  rate 
and  the  velocity,  the  velocity  decreasing  with  decreasing  burn  rate.  Whether 
this  tabulation  Is  some  Indication  of  the  relationship  between  erosive  thres¬ 
hold  and  burn  rate  is  speculation,  but  ought  to  be  kept  In  mind. 

On  the  basis  of  this  revised  interpretation  of  the  CAROE  work,  it  would 
be  concluded  that  triggering  is  made  more  difficult  not  so  much  by  lowering  the 
operating  pressure,  but  rather  by  maintaining  an  adequate  port  flow  velocity 
during  burn.  This  adequate  velocity  may  be  preventing  flow  reversals  or  a  fall 
below  erosive  thresholds. 


2-10 


m 


2.2  EXPERIMENTS  TO  CHARACTERIZE  HETEROGENEITY  EFFECTS 

2.2.1  Prnpellants 

In  attempting  to  characterize  the  heteroaeneity  of  solid  propellants, 
three  propellants  were  especially  formulated  to  highlight  particle  size  effects. 
Each  contained  25%  HTPB  binder  and  751  AP  oxidizer.  The  oxidizer  size  fractions 
were  sieved  to  narrow  each  distrihutinn.  One  propellant  contained  a  monomodal 
fine  (45  ^m)  oxidizer  size,  one  a  monomodal  coarse  (180-212  p.m)  size,  ana  one  had 
a  50/50  bimodal  h'end  of  the  two. 

2.2.2  Microwave  Burner  Tests 


Oynamic  burning  rate  data  with  improved  si gnal -to-noi se  characteristics 
were  obtained  for  the  three  propellants  at  a  sinqle  pressure  of  3.5  MPa  (500 
Psi).  The  reduced  nnise  was  obtained  by  cuttinq  down  on  the  electrical  noise 
superimposed  on  the  recorded  microwave  phase-shift  test  siqnal.  The  recorded 
data  will  be  Fourier-analyzed  to  determine  their  frequency  spectral  density. 

Pressure-coupled  response  function  data  (at  the  same  burning  pressure) 
were  obtained  for  the  three  propellants  at  narrow  frequency  intervals  over  the 
range  40-1500  Hz.  The  dynamic  burninq  rate  frequency  soectra  will  be  compared  to 
the  comhustion  response  function  vs.  frequency  characteristics  to  see  if  some 
dominant  frequency  or  frequencies  appear  in  the  combustion  response.  Microwave 
respnrse  function  data  for  several  other  prnpellants,  obtained  under  an  AFRPL- 
sponsored  proqram,  are  also  available  to  supplement  this  research. 

A  few  additional  dynamic  burninq  rate  and  response  function  measurement 
tests  will  be  made  at  varyinq  mean  hurninq  rates  (varyinq  mean  pressure)  to 
verify  consistency  (higher  mean  burn  rate  should  evoke  higher  frequencies). 

2.2.3  EOAX/SEM  Studies 


This  year's  effort  concentrated  on  the  microwave  dynamic  burninq  rite  and 
response  function  measurements,  and  limited  work  was  done  with  the  EOAX/SEM 
experimental  approach  for  characterizing  the  heterogeneity  of  actual  composite 
propellants,  mainly  in  the  area  of  improvement  of  sample  preparation  techniques. 
The  feasibility  of  this  approach  was  demonstrated  last  year  from  tests  carried 
out  on  the  composite  propellant  formulation  A- 1 3 .  It  was  concluded  that  the 
propellant  formulation  did  fluctuate  on  the  macroscopic  scale  and  these  fluctua¬ 
tions  had  significant  amplitudes  and  considerable  frequency  content. 


In  the  future  EOAX  photographs  will  he  taken  of  the  three  specially  for¬ 
mulated  propellants,  varyinq  maqnification ,  sample  location,  and  sample  orien¬ 
tation.  Frequency  distributions  derived  from  Fourier  analysis  of  the  composi¬ 
tional  fluctuations  will  be  studied  for  ordering  and  for  relationships  tn  char¬ 
acteristic  dimensinns  of  the  fine  structure.  A  statistical  analysis  of  the 
results  will  be  carried  out  to  determine  the  validity  and  consistency  of  the 
appearance  of  preferred  frequencies.  The  relationships  between  any  dominant  fre¬ 
quencies  in  the  propellant  heterogeneity  and  dynamic  burning  rate  data  and  the 
measured  response  function  -  frequency  characteristics  will  be  determined. 
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2.3  HIGH  FREQUENCY  COMBUST  ION  RESPONSE  MOOELING 


2.3.1  Limits  of  the  Quasi -Steady  Gas  Assumption 

The  upoer  frequency  limit  to  the  validity  of  th<=  quasi -steady  qas  as¬ 
sumption  is  a  subject  of  some  confusion.  Numbers  are  often  stated,  but  these 
numbers  are  at  such  variance  that  their  factual  bases  are  suspect.  Therefore, 
as  a  first  step  in  this  work,  it  was  deemed  necessary  to  verify  the  upper  fre¬ 
quency  limit,  and  with  specific  reference  to  the  nitramine/minimum-smoke  pro¬ 
pellant  family  that  is  of  interest. 

The  quasi-steady  qas  assumption  is  valid  if  both  of  the  followinq  condi¬ 
tions  are  satisfied  (Ref.  1C).  First,  the  '•esidence  time  of  the  Qas  in  the 
flame  must  be  small  compared  to  the  characteristic  thermal  time  of  the  solid 
phase.  Second,  the  residence  time  of  the  qas  in  the  flame  must  be  small  com¬ 
pared  to  the  period  of  oscillations. 

The  first  condition  may  be  expressed  as  the  ratio 

b  s  r  X*  cq2  P 

1  *s  cs  RTF  (12) 

Assuming  t he t  r  varies  with  the  0.7  power  of  pressure  and  with  the  pro¬ 
pellant  heat  of  explosion  (Ref.  38),  that  X*  may  be  determined  from  a  first- 
order  Arrhenius  kinetics  law,  and  that  cs  =  cq,  Eo.  (12)  may  be  rewritten 
as 


where  the  subscript  REF  denotes  a  reference  condition,  namely  a  propellant  of 
averaqe  energy  level  burninq  at  6.8  MPa.  Usinq  values  for  the  combustion  con¬ 
stants  that  are  used  in  current  theories  of  burning  of  enerqetic  binder  propel¬ 
lants  (Ref.  39),  a  reference  flame  height  calculated  from  the  theory,  and  tak- 
inq  the  flame  temperature  to  be  that  of  the  primary  or  fizz  zone  (which  is  much 
lower  than  that  of  the  final  flame)  (Ref.  38),  b^  may  be  calculated  to  have 
the  followinq  extreme  values: 

high  energy  propellant  at  a  boost  pressure,  o^  =  0.022 

reference  condition,  =  0,815 

low  energy  propellant  at  a  sustain  pressure,  =  0.011 

Therefore,  the  first  condition  is  always  satisfied.  The  second  condi¬ 
tion  may  be  expressed  as  the  ratio 


(14) 
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h  =  X*  Pf 
2 


which  may  similarly  he  rewritten  as 
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(^ref)  ( 

pre  rj 

(15) 


For  the  specification  that  h2  cannot  exceed  0.1,  the  following  upper- 
limit  frequencies  are  calculated 


high  energy  propellant  at  a  boost  pressure,  f  =  16000  Hz 
reference  condition,  f  -  5100  Hz 

low  energy  propellant  at  a  sustain  pressure,  f  =  900  Hz. 


Thus  the  guasi-steady  gas  assumption  has  a  broader  range  of  applicabili¬ 
ty  than  Is  often  stated,  and  in  many  cases  covers  the  high  freguencies  encoun¬ 
tered  in  motors.  Limiting  the  assumption  to  1000  Hz  is  an  understatement,  par¬ 
ticularly  since  the  phenomenon  is  of  greater  interest  with  higher  energy  pro¬ 
pellants  and  at  higher  pressures. 


2.3.2  Applicability  of  the  Classical  Theory  at  Acceptable  Frequencies 

With  this  guestion  settled,  the  classical  theory  of  Denison  and  Baum 
(Ref.  40)  was  examined  to  see  if  the  numerics  would  be  viable  for  nitramine/ 
minimum  smoke  propellants.  It  was  shown  that  they  are  not  for  A?-composite,  or 
reduced  smoke,  propellants  (Ref.  1).  These  two  classes  of  propellants  differ 
in  the  values  of  their  essential  combustion  constants.  Calculations  of  the 
so-called  "A  and  B  parameters"  in  this  model  were  made.  The  range  of  values  of 
A  (15-20)  and  B  (O.fi-1)  are  such  as  to  generate  satisfactory  response  function 
curves  wih  strong  peaks  and  reasonable  high  freguency  responses,  rather  than 
flat-looking  curves  with  weak  responses.  Therefore,  it  is  considered  that  this 
theory  is  worthy  of  further  us«  in  studying  minimum  smoke  propellants  up  to  the 
freguencies  indicated  above.  This  was  left  for  future  work. 


2.3.3  Reformulation  of  the  Tien  Theory 


There  may  be  difficulties  in  Tien's  work  in  that  solution  of  the  coupled 
second-ordpr  differential  equations  is  not  a  simple  matter  and  often  incurs 
numerical  convergence  prohlems.  Another  aspect  of  the  problem  that  could  bear 
strongly  upon  the  types  of  solutions  obtained  is  his  statement  of  the  boundary 
condition  at  the  flame.  In  this  regard  Tier,  used  a  single  flame,  whereas  it  is 
known  that  the  combustion  of  minimum  smoke  propellants  in  the  gas  phase  is  at 
least  two-staged  (Refs.  41-44),  Therefore,  the  approach  that  was  taken  was  to 
see  how  the  eguations  could  be  simplified  to  avoid  potential  numerical  prob¬ 
lems.  examine  and  perhaps  modify  the  boundary  condition  at  the  flame,  and  in¬ 
corporate  a  two-stage  flame  model.  The  work  would  serve  to  cover  frequencies 
in  excess  of  those  permitted  by  the  quasi-steady  gas  assumption,  and  strive  to 
achieve  more  reasonable  results. 
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Progress  consisted  of  setting  down  the  conservation  equations  for  the 
condensed  phase,  the  primary  flame  zone  and  the  final  flame  zone  in  the  gas 
phase.  The  fizz  zone  was  simplified,  relative  to  what  Tien  did,  hy  assuming  a 
constant  reaction  rate  that  not  only  simplifies  the  energy  eguation  but  also 
eliminates  the  need  for  the  species  eguations.  This  assumption  has  been  justi¬ 
fied  by  the  low  activation  energy  of  the  primary  zone  kinetics  (Ref.  41).  The 
final  flame  zone  is  treated  in  a  simplified  manner  by  assuming  that  a  flame 
sheet  model  applies.  This  assumption  has  been  justified  hy  the  high  activation 
energy  of  the  final  flame  kinetics.  As  a  result,  the  species  eguations  are 
again  not  needed  and  the  energy  equation  is  agair.  simplified.  It  is  considered 
that  the  essence  of  the  problem  is  still  retained. 

Discussions  concerning  the  boundary  conditions  at  the  flames  were  initi¬ 
aled  with  Dr.  Fred  Ctlick  at  the  California  Institute  of  Technology.  This 
aspect  of  the  work  was  not  completed.  Dr.  Culick  expressed  interest  in  this 
work  and  would  be  willing  to  contribute  if  it  is  resumed  in  the  future. 

Further  work  on  this  task  was  delayed  in  favor  of  the  velocity-coupled 
response  function  modeling  work.  This  resulted  from  a  program  reduction,  the 
projected  ability  to  arrive  at  a  certain  level  of  completion  of  the  velocity- 
coupling  work  within  the  known  available  funding,  and  an  awareness  that  Flandro 
(Ref.  24)  was  including  the  high  freguency  regime  in  his  own  theoretical  work 
under  Air  Force  sponsorship.  Depending  upon  priorities,  and  the  nature  of 
Flandro‘s  progress,  it  would  he  proposed  to  resume  this  work  in  the  future. 
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TECHNICAL  PUBLICATIONS 


The  following  publications  appeared  in  the  A I AA  Journal  during  the  past 

year: 

(1)  Sfand,  L.  0.,  and  Cohen,  N.  S.,  "Porous  Plate  Analog  Burner  Study 
of  Composite  Solid  Propellant  Flame  Structure,"  A I AA  Journal,  Vol . 
20,  Mo.  4,  April  1982,  pp.  569-570. 

(2)  Cohen,  N.  S.,  and  Strand,  L.  0.,  "An  Improved  Model  for  the 
Combustion  of  AP  Composite  Propellents,"  AIAA  Journal,  Vol.  20,  No. 
12,  Dec.  1982,  pp.  1739-1746. 

The  following  revised  paper  has  been  accepted  for  AIAA  Journal  publica¬ 
tion,  and  publication  is  pending: 

(3)  Cohen,  N.  S.,  "A  Pocket  Model  for  Aluminum  Agglomeration  in  Compo¬ 
site  Propellants,"  presented  originally  as  AIAA  Paper  81-1585. 

The  following  paper  appeared  in  a  Chemical  Propulsion  Information  Agency 
(CPIA)  publication  during  the  past  fiscal  year: 

(4)  Cohen,  N.  S.,  and  Strand,  L.  0.,  "Effects  of  AP  Size  Distribution 
on  the  Pressure-Coupled  Response  Function,"  1 8th  JANNAF  Combustion 
Meeting,  CPIA  Publication  347,  Vol.  Ill,  Oct.  1981,  pp.  7-20. 
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SECTION  4 


PROFESr  jNAL  PERSONNEL 


The  Principal  Investigator  for  this  program  Is  Leon  D.  Strand  of  the 
Thermochemical  Research  and  Systems  Section  (M/S  125-159),  Jet  Propulsion 
Lahoratory/Cal ifornia  Institute  of  Technology,  48D0  Oak  Grove  Drive,  Pasadena, 
California  91109  (telephone  213-354-3108).  His  co-investigator  Is  Dr.  Norman 
S.  Cohen,  Cohen  Professional  Services,  141  Channing  St.,  Redlands,  California 
92373  (telephone  714-792-88D7 ) . 

Mr.  Strand  has  overall  program  responsibility  and  specific  responsibili¬ 
ty  for  the  experimental  work  performed.  Dr.  Cohen  works  under  a  subcontract, 
and  is  responsible  for  the  analytical  model  developments.  8oth  are  former 
members  of  the  AIAA  Technical  Committe  on  Propellants  and  Combustion  Mr. 
Strand  Is  currently  an  associate  editor  of  the  "AIAA  Journal  of  Spacecraft  and 
Rockets,"  and  Dr.  Cohen  is  a  member  of  the  JANNAF  Combustion  Working  Group. 
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SECTION  5 


INTERACTIONS  (COUPLING  ACTIVITIES) 


5.1  PRESENTATIONS 

The  followinq  presentati cns  were  made  under  this  research  contract: 

(1)  Cohen,  N.  S,,  and  Strand,  L.  0.,  "Non-Steady  Combustion  of 

Composite  Propellants,"  1 982  Joint  AFOSR/AFRPl  Rocket  Propulsion 
Research  Meeting,  Lancaster,  California,  March  19B2. 

(2)  Cohen,  N.  S.,  and  Strand,  L.  n.,  "Combustion  Response  to 

Compositional  Fluctuations,"  A1AA  Paper  83-0476,  21st  AIAA 
Aerospace  Sciences  Meeting,  Reno,  v~\\  ,  Jan.  ig83. 

5.2  INTERCHANGE  FUNCTIONS  WITH  GOVERNMENT  LABORATORIES  AND  CONTRACTORS 

Information  derived  from  past  research  on  the  subject  of  nitramine 
propellant  combustion  was  used  to  complement  current  efforts  in  progress  at 
Thiokol  Corporation  (Dr.  0.  A.  Flanigan,  Huntsville  CM  v . )  under  Air  Force 
sponsorship  and  at  Lockheed  Missiles  and  Space  Co.  (Dr.  G.  A.  Lo,  Palo  Alto 
Research  Laboratory),  The  analytical  model  describing  the  steady-state  burning 
of  solid  propellants  in  general,  developed  in  the  course  of  prior  research,  is 
being  used  to  assist  the  development  of  low  (temperature  sensitivity)  pro¬ 
pellants  at  Thiokol  Corporation  (Dr.  D.  A.  Flanigan,  Huntsvi  1 1  e  Di v. )  unoer  Air 
Force  sponsorship.  This  model  also  has  been  used  in  the  course  of  informal 
contacts  with  other  model  developers  (J.  P.  Renie,  Purdue  University,  Dr.  M. 
W.  Beckstead,  Brighom  Young  University,  M.  K.  King,  Atlantic  Research  Corpora¬ 
tion)  to  compare  certain  model  calculations.  Several  useful  discussions  were 
held  with  Dr.  Beckstead  regarding  the  thermochemical  hases  for  the  energetic 
binder  combustion  model  developed  by  him  under  Air  Force  sponsorship. 

Productive  interchanges  have  also  taken  place  regarding  the  subjects  of 
the  current  research  effort.  R.  R.  Miller  (Hercules  Inc..  ABL  Div.)  has 
observed  and  measured  dynamic  hurning  at  constant  pressure  by  microcinematogra¬ 
phy  of  propellant  strands  containing  coarse  AP.  Additional  evidence  of  multi- 
peaked  response  function  curves  has  been  acquired  at  this  laboratory  and  at  the 
Georgia  Institute  of  Technology  (Dr.  B.  T.  Zinn)  in  related  work  under  Air 
Force  sponsorship.  The  effect  of  AP  size  distribution  on  combustion  response 
has  also  been  the  subject  of  communications  with  Thiokol  Corporation 
(R.  0.  Hessler  and  J.  0.  Hightower,  Huntsville  Div.,  and  D.  P.  Clark,  Wasatch 
Oiv.),  Purdue  Unversity  (Dr.  R.  L.  Glick,  Prof.  J.  R.  Osborn  and  J.  P.  Renie), 
Aerojet  Tactical  Systems  (M.  J.  Oitore)  and  United  Technologies  Corp.  (Dr. 
T.  P.  Rudy,  Chemical  Systems  Oiv.);  two  of  these  involved  practical  motor 
development  programs.  The  EOAX  and  microwave  burner  techniques  may  find  useful 
application  in  the  quality  control  evaluation  of  Space  Shuttle  SRM  propellant 
batches  as  a  result  of  discussions  with  Thiokol  Corporation  (L.  H.  Sayer, 
Wasatch  Oiv.)  and  NASA  (Dr.  B.  Shackelford,  Mashall  Space  Flight  Center).  With 
respect  to  velocity-coupled  instability,  the  reported  work  at  the  Princeton 
Combustion  Research  Laboratory  (Dr.  M.  Ben  Reuven)  has  been  reviewed  and 
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commented  upon,  information  generated  by  M.  K.  King  (Atlantic  Research  Corpora¬ 
tion)  and  K.  Kuo  (Pennsylvania  State  University)  under  Air  Force  sponsorship 
has  been  used  in  this  program,  and  frequent  discussions  have  heen  held  with 
J.  N.  Levine  (Air  Force  Rocket  Propulsion  Laboratory,  AFRPL)  regarding  progress 
under  this  and  related  AFRPL-sponsored  programs.  It  is  planned  to  take  an 
active  role  in  defining  th?  integration  of  the  various  velocity-coupling  ele¬ 
ments  being  worked  upon  in  the  course  of  future  work. 
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AP  ammonium  perchlorate 
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JANNAF  Joint  Army  Navy  NASA  Air  Force 
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SRM  solid  rocket  motor 
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